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Abstract 

First and second order transport coefficients are calculated for the strongly coupled M = 4 SYM 
plasma coupled to massless fundamental matter in the Veneziano limit. The results, including 
among others the value of the bulk viscosity and some relaxation times, are presented at next-to- 
leading order in the flavor contribution. The bulk viscosity is found to saturate Buchel's bound. 
This result is also captured by an effective single-scalar five-dimensional holographic dual in the 
Chamblin-Reall class and it is suggested to hold, in the limit of small deformations, for generic plas- 
mas with gravity duals, whenever the leading conformality breaking effects are driven by marginally 
(ir)relevant operators. This proposal is then extended to other relations for hydrodynamic coeffi- 
cients, which are conjectured to be universal for every non-conformal plasma with a dual Chamblin- 
Reall-like description. Our analysis extends to any strongly coupled gauge theory describing the low 
energy dynamics of N c S> 1 D3-branes at the tip of a generic Calabi-Yau cone. The fundamental 
fields are added by means of 1 ^ Nf -C N c homogeneously smeared D7-branes. 
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1 Introduction 



Hydrodynamic models provide a fairly accurate description of the large-time evolution after 
thermalization of the QCD plasma produced at RHIC [1]. In this respect, a first principle 
computation of the transport coefficients is both relevant and challenging, due to the strongly 
coupled nature of the system. Surprisingly, some features of the known transport coefficients 
appear to be common to many plasmas of strongly coupled theories, including QCD. The 
prototype example is the shear viscosity, whose ratio with the entropy density has, in theories 
with two-derivative gravity duals, a universal value which is quite close to the QCD one [2]. 
This observation justifies a careful analysis of first and second order transport coefficients in 
the class of four dimensional quantum field theories admitting a gravity dual. 

Hydrodynamic studies by means of the gauge/gravity correspondence have been focusing 
to-date on theories without fundamental flavors (apart from what concerns the universal 
shear viscosity [3, 4, 5, 6, 7]). Fundamental flavors are clearly relevant in the RHIC plasma 
[1, 8]. In this paper, the results for the flavor contribution to the bulk viscosity and to some 
relevant second order transport coefficients (relaxation times, k, k*) of the strongly coupled 
M = 4 SYM plasma are presented. 

These results are easily extended to M = 1 plasmas describing the thermal low energy 
dynamics of N c 3> 1 D3-branes at generic Calabi-Yau cones over five-dimensional compact 
Sasaki-Einstein manifolds X 5 (the Af = 4 SYM case corresponds to X 5 = S* 5 ). These are 
conformal plasmas without fundamental fields. Flavors are added by means of Nf 3> 1 ho- 
mogeneously smeared D7-branes extended along the radial direction and wrapping compact 
three- manifolds X 3 in the space transverse to the D3-branes [9, 5, 10, 7]. We will just focus on 
the case in which all the flavors are massless; the related non-Abelian flavor symmetry group 
is explicitly broken into a product of Abelian factors, due to the smearing. The addition of 
massless flavors induces a breaking of conformal invariance at the quantum level. 

The reason to analyze these theories stems from the fact that, despite being plagued by 
a Landau pole in the UV, they are the simplest examples where thermal flavor effects at 
strong coupling can be reliably studied, providing the first manageable toy models for the 
QCD plasma in its near conformal regime. 

Our results are obtained at next-to-leading order in a perturbative expansion in €h ~ 
X h Nf/N c where is the 't Hooft coupling at the energy scale fixed by the plasma tempera- 
ture. The precise coefficients defining th depend on the model and hence on the volumes of 
the X 5 and X 3 spaces. For the Af = 4 case, for example 



Within this perturbative approach, the flavors can be considered as "deformations" of the 
conformal plasmas. At zero temperature, the flavor superpotential term, which drives the 
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breaking of conformality, can be treated as a marginally irrelevant "deformation" [10], ac- 
cordingly. 

At second order in the gravity solutions derived in [7] provide a completely reliable 
description of the dual field theories in the planar limit and at strong coupling, by accounting 
for the backreaction of the branes supporting the flavor degrees of freedom. Thus, they allow 
for the calculation of the transport coefficients by means of standard procedures. 

The main results are collected in the following subsection. A review of the relevant back- 
grounds can be found in section 2, the main steps of the calculations are reported in section 3 
and further details are given in appendix A. In section 4 we will provide an alternative simple 
calculation for the bulk viscosity, using the arguments presented in [11]. In particular we will 
show how an effective single-scalar five-dimensional holographic model in the Chamblin-Reall 
class [12] captures the leading conformality breaking effects due to the marginally irrelevant 
flavor "deformations" . We will also show that the same approach can be successfully applied 
to cascading plasmas - where conformality breaking is driven by marginally relevant oper- 
ators - at leading order in the perturbative expansion introduced in [13]. We will suggest, 
in turn, that in the limit of small "deformations", there are certain universal relations for 
transport coefficients of gauge theory plasmas (with gravity duals) where the leading con- 
formality breaking effects are driven by marginally (ir)relevant operators. Our proposal is 
a natural extension of the results in [14], valid for relevant or exactly marginal deforma- 
tions. Finally, we will suggest possible universal relations involving the bulk viscosity and 
the interaction measure. As we report in the following subsection, the results of section 4 
allow us to propose a class of universal relations for every non-conformal plasma with a dual 
Chamblin-Reall (effective) description. 

1.1 Main results 

Up to second order, the hydrodynamic expansion of a non-conformal plasma is known to 
be determined by two first order transport coefficients, i.e. the shear and bulk viscosities, 
and thirteen second order transport coefficients [15], the most important ones being the 
relaxation times, which are relevant for numerical hydrodynamic simulations. We refer to 
[15, 16] for the notation of the transport coefficients. 

The bulk viscosity and a combination of the "shear" and "bulk" relaxation times r n , Tn can 
be derived from the dispersion relation of the scalar hydrodynamic modes (sound channel) 



In this equation, u is the frequency of the mode and q its momentum; r\ and ( are respectively 
the shear and bulk viscosities; s, T and c s represent the entropy density, temperature and 
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speed of sound of the plasma. Finally, r e ^ is an "effective relaxation time" 1 which for 
non-conformal plasmas is the combination 

i*" = 4 { ■ (1.3) 
1 + | 

The second order coefficient k and a combination of r n and k* can be derived from the 
retarded correlator of the tensorial mode [17, 15] 

GT V =p-i^+ (rrr, - \ + ^ ~ ^ + 0{q\ c 3 ) , (1.4) 

where p is the pressure. 

In [7] the following quantities were obtained 

r = T„(l-^-^), P = Po(l-^), £ -3p = |e L 

Cs = A=(l-—e(), - = — , (1-5) 

v/3 V 12 V ' s 4vr ' v ; 

where T and po = 8V oilx° 5 ) are ^ ne temperature and pressure of the unflavored conformal 
plasmas, and e is the energy density. 

The new results in this paper, calculated up to C(e^), are 

\ = 2 ^ 

r"'T = r,,„T + ^^4, (1-7) 

— K = -2-Ko, (1.8) 

— k' + jjtv = -^oT^o + ^o — ^ e 2 , 1.9 
p po po \ 8 1287rT / 

2-log2 T 2 1 Tom 1 n in x 

t„, T = — , — K o = ^, = -, (1.10) 

are the corresponding values in the conformal plasmas [17, 18]. 2 

Let us comment these results in turn. The bulk viscosity (up to 0(e 2 h )) saturates the bound 
proposed by Buchel in [19], i.e. (/rj > 2(1/3 — c 2 ). In section 4 we obtain again this result 
in the framework of [11, 14], showing that it is just the source of the marginally irrelevant 



where 



1 This terminology is borrowed from [16]. 

2 In the conformal case t„- = t^q = r e ^ , since ( = 0. 
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operator driving the breaking of conformal invariance (essentially the flavor superpotential 
term dual to the dilaton field) the one which contributes to the bulk viscosity at leading 
order. We see no reasons for the bound to be saturated at all orders. 

The numerical value of (/rj, of order 10~ 3 for ~ (9(10 _1 ) [7], is quite small as compared 
to the expected one in the near-conformal region of QCD, i.e. (/rj ~ O(10~ r ) [20]. This 
is expected and due to the fact that the models at hand compute just the quantum flavor 
effects and do not include the pure YM contribution to the trace anomaly, which is the main 
source to (/rj in QCD. 

The result for k is usually given in the conformal case in the combination T k /so = 
TqK /Apq. In the present case it reads 



Tk, _ T /t 
s s 



(l-f). (1.11) 



As such, all the results (1.6)- (1.9) tell us that the leading corrections to the conformal values 
of dimensionless combinations of transport coefficients consist of a (coefficient dependent) 
constant times a common function of the temperature, 3 which in the present case is e 2 ~ 
log -2 T. This is in agreement with an extension of the results presented in [14] to cases with 
marginally (ir)relevant deformations. In view of the general behavior discussed in section 
4, it would not be surprising if the numerical coefficients in (1.6)-(1.9), written in terms of, 
say, 5 S = 1 — 3c 2 using eq. (1.5), turned out to be universal in all the marginally (ir)relevant 
deformations of 4d conformal theories with two- derivative gravity duals. 

Unfortunately, T n , Tu and k* cannot be disentangled in the computations reported in this 
paper. Nevertheless, the results in (1.7), (1.9) give very interesting indications. Comparing 
(1.3) with (1.7) and (1.9) using the notation 4 

Ttt = t Wj0 + T^^eh + r n>2 el, r n = r n ,o, k* = «S + k\ e h + tC 2 e 2 h , (1-12) 

where in each of these expressions it is understood that we neglect higher order terms in eh, 
it follows that 

1 r no _ 17 16-tt 2 

r.,i-g^o, ^,2 + ^-- 128^.0+128^ ' { ' 

* * n L ^2 19 16 + tt 2 

Ko = Kl = °' T - 2 + ^ = 3^ r -°-T28^T-- (L14) 

In section 4 we will note that some theories where conformal invariance is broken by 
marginally (ir)relevant operators are effectively described, for what concerns the bulk vis- 
cosity and speed of sound and at leading order in the deformation, by dual Chamblin-Reall 



3 We thank A. Cherman, T. Cohen and A. Nellore for this observation. 

4 Sincc Q/rj ~ ef, only the zcro-th order term in m is relevant in our perturbative analysis. 



5 



models, i.e. metric plus single scalar theories with exponential potential. For this class of 
models, in [21] it was found a relation between r n , rn which reproduces exactly the result in 
(1.13). 5 This means that the theories analyzed in this paper, where conformal invariance is 
broken by marginally (ir)relevant operators, are effectively described, at least at third order 
in the sound channel and at leading order in the deformation, by Chamblin-Reall models. 

Crucially, the latter include also the non-conformal theory in [22], where the conformal 
transport coefficients determine all the others, in particular rn = r n and k* = — ^(l — 3c^) 
[15]. For the D3-D7 plasmas analyzed in this paper, postulating rn = r n (i.e. rn,o = Tn,o) 
implies from (1.13), (1.14) that k\ = —/to/4, that is, using the value of the speed of sound 
in (1.5), precisely the relation above. 

Thus, ignoring the possibility of a mere coincidence, we are led to conjecture that for 
all the theories effectively described by a Chamblin-Reall model, the transport coefficients 
satisfy the relations reported in [15]: 

A 4 = 0, 

(1.15) 

6 = y(l-3a, £4 = 0. 

In turn, the relations (1.15) allow to make a prediction for all the second order transport 
coefficients, apart from Ai, A2, A3, up to 0(e\) for the D3-D7 plasmas (using the same notation 
as in (1.12)): 

Kq — K 1 — U , 

r n,o = Ttt,o , 
£2,0 = £2,1 = , 

The results in section 4 suggest that similar predictions can be made for the cascading 
plasmas. Note that r n = r w would imply that r e ^ = r n . Moreover, from (1.7) it would 
follow that r w T > TnfiTo at order 0(e 2 h ). So, these results would support the conjecture 
in [16] that T ni oT is a lower value for the relaxation times in theories with two-derivative 6 
gravity duals. 7 One could also conjecture, analogously, that kq constitues an upper bound 
for the same class of theories. 

5 We thank T. Springer for pointing out this agreement. 

6 Beyond two derivatives, as for the shear viscosity, while finite coupling corrections enhance the relaxation 
time [23] , finite iV c corrections reduce it [24] . 

7 Similar bounds on the speed of sound were proposed in [25], [26], and a possible universal relation among 
some second order coefficients were proposed in [27], [28]. 
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* K * * n * '7T,0 -v n 

^2 = -^p ^,0 = ^,1 =0 ' r -,2 = g"> A 4 = 0, 

c Vor^fi _ _ _ A 3i0 t _ n 

?2,2 — 2 „ , ?3,0 - ?3,1 - U , ^ 3 ' 2_ ^g" _U ' ?4 — U - 



Finally, we note that the known results [17, 18, 22, 21], including the ones of this paper, 
are compatible with an analog of Buchel's bound of the form 




(1.17) 



where d is the number of spatial directions of the plasma (d — 3 for the models considered 
in this paper). Again, analogous bounds could be proposed from the relations (1.15). It is 
necessary to underline that these proposals are not based on first-principle arguments and 
they would clearly require further, more solid, confirmations. 

In conclusion, the flavor effects are found to enhance both the bulk viscosity (1.6) (simply 
because, breaking conformal invariance, they produce a non-zero () and the "effective relax- 
ation time" (1.7), while they reduce k in the usual combination (1.11) but do not modify 
it in the combination (1.8). Finally, breaking conformality, the fundamental flavors should 
give non-trivial coefficient r n , k* (1.13), (1.14). 8 Unfortunately their values cannot be reli- 
ably disentangled from r w with the present gravity analysis, but we believe they satisfy the 
relations (1.15), which would imply that T n is enhanced by flavor effects. 



2 The background 

In this section we summarize the results obtained in [7] for the non-extremal gravity solution 
describing an intersection of two sets of N c 1 D3-branes and 1 <C Nf <C N c D7-branes, 
where the backreaction of both stacks is taken into account. When no D7-branes are present 
the solution is AdS 5 x X 5 with a black hole. The presence of backreacting D7-branes intro- 
duces a squashing in the internal X 5 (described as a U(l) fibration over a four-dimensional 
Kahler-Einstein (KE) base) and the Einstein frame metric reads 9 

ds 2 = h~ l/2 [-bdt 2 + dxidxi] + h 1/2 [S^F^dr 2 + r 2 (S 2 ds 2 KE + F 2 (dr + A KE ) 2 )] . (2.1) 

Here dA KE /2 = J KE is the Kahler form of the four-dimensional base of X 5 . Moreover 

r 4 ' r 4 ' 4 Wol(X 5 ) ' V ; 

where the subscript h stands for "horizon" and 

p = ! e *fi| 2r4 ~ r ^ | e * ( 17 94 2^-^ 5 (2^-^ 



24 V. 6r 4 -3r 4 / 1152 V. 9 2r 4 - r 4 9 (2r 4 - r 4 ) 2 
r 8 h (rf-r A ) , Q1 r 



9 (2ri 



„4\3 



- 48 log • 



8 At least one of them is non-zero. 

9 With respect to [7]'s notation we have removed the tildes from F, S. 
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up to second order in the perturbative expansion parameter e* = e(r*), where 



Q f e* 



Vol{X 3 



(2.3) 



(2.4) 



16ttV ol(X 5 ) N c 

and A = Ang s N c e®; here ^o/(X 3 ) is the volume of the three-dimensional compact mani- 
folds wrapped by the homogeneously smeared D7-branes. For the M = 4 plasma we have 
Vol(X 5 ) = 7T 3 , Vol(X 3 ) = 2tt 2 . 

The D7-brane sources induce a running of the dilaton $(r), which reads 



+ e* lo, 

9 



+ - 



Li- 



r K 

h — 

n 72 



r 



1 - 7T4 \ + 121 °g- + 361 °g 2 - + 

2rt - rt r* 



(2.5) 



The non-trivial RR field strengths on the background are 

F {5) = Qc (1 + *)e(X 5 ) , F ( i) = Q, (dr + A xs ) dF {l) = 2Q f J KE , (2.6) 

where e(X 5 ) is the volume element of the internal space. 

As mentioned above, at order zero in e* the solution is the standard AdS$ x X 5 black 
hole. The UV cutoff A*, corresponding to the radial position r*, is ensured to be well below 
the Landau pole Alp (mapped, in turn, to the position r = r^p at which the exact solution 
for $(r) blows up) if e* <C 1, and thus represents the scale up to which the solution is 
under control. At A* a UV completion is needed. We focus here on the IR properties of 
the system, such that we can safely discard contributions coming from the UV completion, 
which are power-like terms suppressed as (A/p/A*)™. Thus, in the computations below we 
systematically discard the (subleading) corrections in r h jr\. In this regime the previous 
expressions for F, S and $ get the simpler form 



F = 1 



eh 



17 



1 + 



24 + 1152 

z_h 1 
24 128 

r 



— log — 

24 ^r h 



F h 



— log — 

24 & r h 



€ h i„„ r — c , € h 



el . r el . , r e? 



$ h + e h log — + log — + -£ log 2 — + -^Li 2 
r/j 6 r h 2 r h 16 



r 4 



Notice that, being the temperature T proportional to r ft at leading order [7], one has 



.2, r h 

e h = e if + e, log — , 



(ie/j _ 2 
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(2.7) 



(2.8) 



8 



which implies that while e runs with the energy scale (and thus with the temperature), e 2 is 
constant if we neglect e 3 and higher order terms. The interaction measure (e — 3p)/T A and 
the shifted speed of sound c 2 — 1/3 (see eq. (1-5)), which are zero for conformal plasmas, 
scale like e 2 as an effect of quantum conformality breaking due to the dynamical massless 
flavors [7]. Those observables are thus independent on the temperature at next-to-leading 
order in e h . We will come back to this observation in section 4. 

2.1 Reduced action 

In the following we will study hydrodynamic transport properties of the plasmas described 
above. We first rewrite our ten-dimensional model as an effective five-dimensional one by 
integrating over the internal manifold X 5 . The five-dimensional reduction of the action reads 
[10] 



_ Vol{X b ) f 5 r ~— 

^io J 



R[g] - f(df) 2 - 20(dw) 2 - l -m 2 - v(®, /, w) 



, (2.9) 



where g mn is the five- dimensional metric, R[g] is its scalar of curvature and the potential 
describing the interactions of the three scalars reads 

V(<f>, f, w) = 4ef f+2w (e 10w - 6 + Q f e") + l -Q)e^^ w+2 " + ^ef f . (2.10) 

Hereafter we set R — 1 and a' = 1 for simplicity. In these units Q c = 4 as can be read from 
(2.2). 

To get the above results one starts from a reduction ansatz of the form 

d s \ Q = effg mn dx m dx n + e-^ f+w) ds 2 KE + e 2 ^ w ~ f \dT + A KE f . (2.11) 
On the background (2.1) we have 

/ = -1 log (S 4 F) , w = 1 log (F/S) . (2.12) 



The five- dimensional non-extremal background metric reads 



ds\ = r 2 e ^ f [-bdt 2 + dxidxi] + e"^ 1 ^ 1 ^ = -c^dt 2 + c 2 x dxidxi + c 2 R dr 2 . (2.13) 



From the dimensional reduction we see that we have three scalar fields. In the perturbative 
expansion in e* our models (both in the extremal and in the non-extremal case) can be seen 
as "deformations" of the unflavored AdS$ x X5 (BH) solutions. The AdS background is a 
minimum of the potential at order zero in e*, in which case / = and w — 0. The field 
/ is dual to an irrelevant operator of dimension A = 8 whose form is TrF 4 . It drives a 
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deformation from the AdS 5 x X 5 to the non-near horizon D3-brane background. The field w 
is dual to a vev for an irrelevant operator of dimension A = 6. It is of the form Tr(WQ,W Q ) 2 
and is the responsible of the squashing of the transverse four-dimensional Kahler-Einstein 
base and the fibration. The dilaton $ is dual to the insertion of a marginally irrelevant 
operator, actually the flavor term in the field theory (T = 0) superpotential [10]. This is the 
source term which is responsible for the breaking of conformal invariance at the quantum 
level. 

Notice that the action (2.9), and so the equations of motion for some of the perturbations, 
coincides with the one in [29] (but for the definition of the potential). The equations of motion 
for the fluctuations not coupling to the scalar fields (that is, insensitive to the potential) can 
be read in [30]. 



3 Fluctuations 

Following the standard procedure [31], we will consider fluctuations of the fields present in 
the five-dimensional action. It is easily seen from the expansion of the DBI action that 
the perturbations considered in the following do not couple to those of the flavor branes 
at the linearized level, thus the relevant set of fluctuations we need to consider is ^(r) — > 
tt(r) + 5t>(x», r), with = {g mn , $, /, w}. 

We will assume that the perturbations take a planar wave form in Minkowski space, 
<5#(^,r) = e-^-^ir). Thus, these fields can be classified according to their trans- 
formation under the little group 5*0(2), which is the remaining symmetry of the system 
(rotations in the x — y plane). We define 

Sg tt (x^r) = e-^-^cU^Huir), (3.1) 

Sgmnix^r) = e-^-'V^WTUW, (m,n)^(t,t) (3.2) 

5<f>(x»,r) = e-^ 9Z V(0, (3.3) 

6f(x",r) = e- i(uJt ' qz) B(r), (3.4) 

6w(x",r) = e- i{ut - qz) C{r) , (3.5) 

choosing the gauge Ti. rm {r) = 0. The classification of the different channels gives a tensorial 
mode (7i xy ), vectorial modes (Titx, 7~C ZX , Tity, 7~Lzy) and scalar modes (TCu, 7~t±± = Ti. xx + H y y, 
7~Lzz, "Htzi ¥■> C)- Each kind of perturbation can be expressed in term of gauge invariant 
quantities under the residual gauge symmetry [31, 29, 30] 

Tensorial — * Zt = Ti, xy , 
Vectorial — > Zy = qHt x + ul~L zx , 



Scalar -> Z s = 2H ZZ + 4-H tz 

LO 



q 2 c' t Ct 
u 2 c' x c x 



q cfi 

T~t±± + 2— -o-Htt , 

u 2 c 2 x 
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!og [4] 
/' 

B — ; r , t H±± 



log 



w 

log [41 

In all the cases, studying the differential equations at the horizon we find that the solutions 
behave as — l) ±l "25^, where to = o;/(2r^) and = 7rT . We choose the index with 
negative sign to have incoming boundary conditions at the black hole horizon. Moreover, as 
in [32], we impose that the fluctuations vanish at the UV cutoff scale related to r*: crucially, 
our results will turn out to be independent of r* up to suppressed terms in powers of r\/r*. 

In the following we are going to study the tensorial and scalar perturbations, while the 
study of the vectorial perturbation is relegated to appendix A. 



3.1 Tensorial perturbation 

We scale vo — > \hyd&>, 9 —> ^h y d% where q = g/(2r/ l ) and \h y d is a parameter keeping track 
of the order of the hydrodynamic expansion. Define 

4 ^ 2 2 

^ = Ct(i-^) ££4'*^, (3.6) 

where higher order terms in e* and \h y d, which we will not study, are not taken into account. 
The equation satisfied by the perturbation is reported in appendix A. As expected, Dirichlet 
conditions cannot be imposed, showing the absence of a dispersion relation in this channel 
[31]. However one can write the hydrodynamic expansion of the retarded correlator. To 
do this we have to evaluate the action on-shell. This action is singular when evaluated 
at r = r* — > oo (it goes as r*). To cure this divergence we have to add the following 
counterterms [33] 

v,,,,, ^, ul k + ^t^ f'^^ — K • ftt^ — fWo V o /,' - ) , (3.7) 

IK W J ^ K 1Q J ■ Z ' 



where K is the scalar associated to the extrinsic curvature, 7 is the four-dimensional metric 
at the boundary, C[4>] is a function of the scalars and W[0] is the superpotential 

W[f, w, $] = e§ / [Q c e 5/ + Q f e f ~ Aw+ * - Ae f+6w - Qe f - 4w ] , (3.8) 

from which the potential (2.10) can be derived as 



2 



3 fdw \ 1 1 (dwV fow 



80 V df J 40 V dw J \d$ 



- \W 2 . (3.9) 
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The function C[4>] satisfies the differential equation [33] 



1 1 

2 ~ 4 



3 dWdC 1 dWdC dW dC 
+ T77 t; + 



80 9/ 9/ 40 <9w «9u> <9$ <9$ 



+ J-CW = 0. (3.10) 



Although we do not know the exact form of C[(f>], we can extract physical results, since 
the divergence balanced by this function goes in the UV as r 2 , being the next-to-leading 
order suppressed as r\ 2 , i.e., it does not affect the finite part from which the hydrodynamic 
transport coefficients are obtained. For completeness we give its leading behavior, needed to 
cancel the divergence 

0[/,»,*]«l + ^«.-^ + 0(r-). (3.11) 

The Fourier transformed, quadratic-in-fluctuations, on-shell boundary action is 

S = V °^ ) f d A kH_ k J^{k,u)H k) (3.12) 

with H k the boundary value of the fluctuation. The retarded correlator of the corresponding 
components of the energy momentum tensor 

G x ^ xy (uj, q) = ~ij dtd^e^-^emT^x),^^)]) , (3.13) 

is related to the on-shell action by [34] 

C%™(u,q) = -21m[F(k,r m )]. (3.14) 

Plugging the solution of the equation of motion for Zt, equation (A.l), in the finite action 
(3.7), it is straightforward to derive the flux JF(A;, r*), and so the correlator G x ^ xy {oj, q). 
Using also (2.8), we get 



24 + 19m> - 4q 2 + 2t» 2 (2 + 3tt 2 + 22 log 2) 
192 



4) , (3-15) 



which, compared to (1.4) and using the expression for the temperature in (1.5), confirms 
the values of the pressure p and shear viscosity rj from [7] and gives the new results in (1.8), 
(1.9). 
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3.2 Scalar perturbations 

Focusing now on the scalar fluctuations we write for each perturbabion Za=s,b,c,<p the ansatz 

2 2 2 

w = ^c s , fe e*q-2i^ 7fe e*q 2 + 4^e*q 3 , (3.16) 

k=0 k=0 k=0 



t»T „ „ 
4\^ i ^ i 2 2 



^ = ^(l-i) ' EE^V^, (3.17) 

j=0 fc=0 



where 7fc,£fc are the order coefficients of the adimensional combinations 

7 = 7rTo r, t=(vr T ) 2 I( c ^-^). (3.18) 

Higher order terms in e* and q, which we will not study, are not considered. 

The relevant equations for the perturbations are reported in appendix A. The calculation 
can be performed imposing regularity at the horizon. Once this is obtained, one can ask 
for Dirichlet conditions at the UV cutoff, r*, eventually taking the limit r* — > oo. Here 
we present only the results relevant for the physical observables. The solution is given in 
appendix A. With the Dirichlet condition at the boundary we find 

1 1 

Cs '° ~~ v^' c s,i-°> °s,2 - ~ 12v /3 ' {S.19) 

1 1 17-161og[£] 

7o = 77 , 7i = 77j > 72 = 



6 ' ' 48 ' ' 768 

3-21og2 3-21og2 , 57 - 3vr 2 - 22 log2 - 24(3 - 2 log 2) log[£] 



^0 — ~ — 7= i ti — ^= , ti — 



24v/3 ' 96^ ' 2304^ 

which confirms the result for the speed of sound found with the thermodynamics in [7] and, 
using (1.2), (1.5) and (2.8), gives the new results reported in (1.6), (1.7). 



4 A simple perturbative approach to first order hydro- 
dynamics 

Let us consider a simple five- dimensional gravity model with a minimally coupled scalar field 





5,5 = Y 1 I d5x ^/~ det 9 

5 J 



R[g]-^d<t>) 2 -V{<t>) 



(4.1) 



and assume that this model admits an AdS (black hole) vacuum - dual to a four- dimensional 
(thermal) CFT - when the scalar field is turned off. Considering thermal cases in which <fi is 
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dual to a source for a relevant operator of asymptotic dimension 2 < A < 4, the authors of 
[11] provided a simple expression relating the bulk viscosity 10 of the dual field theory plasma 
with the five-dimensional scalar potential: 



C - \hV>(6 )\ 2 ( V ' {<Ph) 



(4.2) 



An analogous approximate expression for c 2 — 1/3 was proposed in [35]. In (4.2) the prime 
means derivative w.r.t. and <p h is the value of the field at the horizon. The coefficient 
h-ni^h) is determined by solving the equation of motion for the SO (3) invariant fluctuation 
Hn = H22 = i?33 = e lujt hn((f)) of the three-dimensional spatial metric components 11 at 
uj = 0. In the r = gauge, in which the five- dimensional radial coordinate is identified with 
the background scalar field 0, the equation is given by 



v 



6A> 



-Y' 



h 



11 



(4.3) 



where the functions A, Y enter the background metric as 



ds 2 = e 2A [-bdt 2 + dxidxi] + e 



2Y' 



(4.4) 



In order to solve eq. (4.3) one imposes regularity at the horizon (selecting only incoming 
waves in the uj 7^ case) and the (normalized to fix the residual three-dimensional scale 
factor) boundary condition hu — > 1 at the asymptotic AdS boundary. 

In [14] (see also [36]) it was shown that the coefficient hf}(<f) h ) depends only on A (and not 
on the details of the five-dimensional potential) at leading order in a perturbative expansion 
around the conformal background (i.e. in a large T limit if the deformation related to is 
relevant). The general results of [14, 36] can be eventually extrapolated to exactly marginal 
deformations, in which case they consistently give trivial hydrodynamic coefficients (i.e. 
C = c 2 — 1/3 = 0); however, they do not apply to theories where conformality breaking is 
driven by marginally (ir)relevant operators, like, notably, Yang-Mills as well as the D3-D7 
and the cascading plasmas [13]. In [37] it was shown that, at T > T c , hf}{(f) h ) — > 1 in a 
five-dimensional phenomenological holographic dual to Yang-Mills. In the following we will 
show that the same result applies to the D3-D7 and the cascading plasmas in analogous 
asymptotic regimes. 

It is relevant to notice that, as it was shown in [11], eq. (4.3) always admits a constant 

1 to satisfy the boundary conditions) for 



solution (which in turns implies that h u (4>h) 

10 This is the only non trivial first order hydrodynamic coefficient, since rj/s = l/(4n) for any strongly 
coupled plasma with a two-derivative gravity dual [3]. 

n In fact the bulk viscosity is related, via Kubo formulas, precisely to the (low frequency limit of the) 
50(3) invariant retarded Green's function of the operator Tn + T22 + T33, where is the field theory 
stress-energy tensor. 
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Chamblin-Reall backgrounds [12], where V — V e 7< ^, with V < and 7 constants. As we will 
see in the following, the 5d duals to the D3-D7 plasmas effectively behave like the Chamblin- 
Reall models, at leading order in the e-expansion around the conformal fixed point. As we 
will show, the same applies to the cascading plasma [13] in an analogous "perturbative" 
regime. It is worth underlining that other relevant non-conformal plasmas [30, 22] have a 
dual gravity description precisely given by a model in the Chamblin-Reall class [35, 11]. We 
argue that this is the reason why the above mentioned systems (at leading "perturbative" 
order, in the case of the D3-D7 and the cascading plasmas) have so many common features 
(for example they all saturate the bulk viscosity bound proposed in [19]) despite having 
different microscopical content. 



4.1 Bulk viscosity of D3-D7 plasmas 

Reducing the D3-D7 models to five dimensions (see eqs. (2.9), (2.10)) we have seen that 
there are three scalar fields in the action. However, only one of them, namely the dilaton, is 
dual to the source for the (marginally irrelevant) deformation driving our theories away from 
conformality [10]. It is thus conceivable that the bulk viscosity, which is turned on when 
conformality is broken, is primarily determined by the dilaton field in an expansion around 
the AdS-BH solutions. 

To better understand the role played by the various scalars, let us consider the quantity 

in three different cases, where we identify <fi with one of the three fields entering in the 
potential (2.10), taking the other two scalars fixed to their background values. At order e\ 
we get 

V />tt = 0, V* = |. (4.6) 

This indicates that the bulk viscosity (and the speed of sound) can be determined in our 
models by considering just the dilaton $ as the "active" field in the game. The other two 
scalars, / and w, do not contribute at leading order and they can be fixed to their background 
values. Due to this observation we can immediately apply the recipes of [11], based on a 
single-scalar five- dimensional model, to our cases. 

Let us define, at first order in e^, a new radial variable 

= $ - $ fe = e h log — =>• — = e£ , (4.7) 

rh r h 

from which we can re-express our five-dimensional metric (2.13), in the form (4.4) with 

A = — + ^(l + 0) + const, Y = ^(1 + 0) + const. (4.8) 
€h 24 6 
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From these expressions we find that the term in h u in (4.3) vanishes (up to 0(e|) terms). 
Just as for the Chamblin-Reall models [11], hu = constant is thus a solution to (4.3) at 
leading order. This implies that h^i(4>h) = 1 to satisfy the boundary conditions. As a result, 
Buchel's bound on the bulk viscosity [19] is saturated. In fact we obtain 



V V{cj> h y 9 



(4.9) 



Notice in turn that c 2 s — 1/3 = — e^/18 (as computed in [7] and in the previous section) can 
be expressed as 

just as it happens in the Chamblin-Reall cases [35]. 



4.2 The bulk viscosity of the cascading plasma 

Many relevant five-dimensional gravity models, dual to well studied non-conformal field 
theories, contain more than one scalar. A well known example is given by the cascading 
conifold theory [38, 39] whose hydrodynamics has been studied in great detail by Buchel 
and collaborators (see for example [40]). This theory describes the low energy dynamics 
of N regular and M fractional D3-branes on the conifold. When M — 0, the theory has 
gauge group SU(N) x SU(N) and it is conformal with an AdS 5 x T 1,1 dual [41]. The 
addition of fractional branes modifies the gauge group to SU (N + M) x SU (N) and breaks 
conformal invariance. The combination g{ 2 — g 2 2 of the gauge couplings, in fact, acquires 
a logarithmic running with the scale g^ 2 — g 2 2 ~ Mlog(fj,/Am). The marginally relevant 
operator TrF 2 — TrF 2 is mapped to a massless scalar field in the dual five-dimensional 
gravity description (see e.g. [42] for a complete scalar/operator map in the deformed conifold 
theory). This is actually the supergravity modulus arising from the integral of Bi over the 
two-cycle of the conifold. The other massless field, the dilaton, is dual to TrF 2 + TrF 2 and 
(differently from the D3-D7 cases examined above) is just a constant on the background at 
T = 0. The five- dimensional effective action reads 

^ = Vd } T 2 1 ' l) [ d 5 x^-detg [R[g] - C km - V(f, w, $, K)] , 

Ckin = f(df) 2 + 20(dw) 2 + i(9$) 2 + -L e -* +4 f +4w (dKf , 

V(f, w, $, K) = Aeff +2w {e Ww - 6) + P 2 e* + f f~ 4w + -K 2 ef * . (4.11) 

In a' — 1 units, P ~ g s M and K is proportional to the effective number of regular D3- 
branes, which is running with the scale if M ^ [38]. At P = the previous action has an 
AdS (BH) vacuum where (setting the AdS radius to one) K — K* — 4. On this background 
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$ = / = w — 0. The fields /, w are mapped to irrelevant operators with A = 8, 6 just as in 
the D3-D7 setup. 

A regular non-extremal solution of the equations of motion following from (4.11), has been 
found using a perturbative approach [13] analogous to the one we have adopted in the D3-D7 
case, namely by means of an expansion in 5 = P 2 /K*. This parameter has a non zero beta 
function at 0(5 2 ): its logarithmic running with the temperature (from which it follows that 
5 <C 1 at large T 3> A IR [13]) can be thus neglected at 0{8). At this order and in our 
units, 5 = P 2 /A. If we rewrite K = 4 + \/2Px we can see that the field x has canonically 
normalized kinetic term and it is massless, around the AdS (BH) vacuum. This is precisely 
the field dual to the TrF 2 — TrF 2 operator mentioned above. 

In the non-extremal case, the functions f,w,Q (resp. x) receive the first corrections at 
0{5) (resp. 0(y/5)). If we now consider the quantity V<j, defined in (4.5) we get 

V />w> » = 0, V x = |f, (4.12) 

at leading order. Precisely as in the D3-D7 case, the only "active" scalar, for what con- 
cerns the leading perturbative contribution to the bulk viscosity, is the one related to the 
marginal (in this case marginally relevant) deformation which sources the breaking of con- 
formal invariance. It is possible to verify, as above, that h^((f>h) = 1 at leading order, such 
that 

C VW 8 

rw = s J ' (4 ' 13) 

nicely reproducing the results found in [40] by means of the alternative analysis of the 
hydrodynamical pole in the stress-energy tensor. Notice, again, that the speed of sound 
c 2 — 1/3 = —45/9 (as computed in [40]) satisfies the relation 

i _ivw (4 14) 

* 3 2 V(4» h )'' ' [ ' 

just as it happens for the Chamblin-Reall models. From the above expressions we see that 
Buchel's bound is saturated at leading order, consistently with the results found in [40]. 



4.3 Comments 

From the above results it is tempting to propose that 



C _ v'^ h ) 2 

V V^h) 2 




(4.15) 



in the vicinity of a conformal fixed point - i.e. at leading order in a "small" (resp. "large") 
T "perturbative" expansion - for every four-dimensional plasma with a five-dimensional 
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two- derivative gravity dual where conformality breaking is driven at the quantum level by 
a marginally irrelevant (resp. marginally relevant) operator. A common feature of this 
kind of plasmas is that the perturbative expansion parameter 5 runs as the inverse of the 
logarithm of the temperature and hence (3t[S\ = T(d5/dT) ~ ±5 2 . This in turn implies 
that (5t[{s — 3p)/T 4 ] is subleading w.r.t. (e — 3p)/T A : the latter (as well as c 2 — 1/3) is thus 
effectively constant at leading order. 

The relations (4.15), easily extended to d ^ 3 space dimensions (just replacing 1/3 with 
\jd on the r.h.s), are precisely satisfied by other known non-conformal models [30, 22] whose 
dual gravity description is in the Chamblin-Reall class. The results of [14, 36], instead, show 
that gravity duals of relevant deformations do not belong to this class, as eqns. (4.15) do 
not hold. 

A well studied example of a plasma where conformality breaking is driven by relevant 
deformations is the M = 2* one [43]. In this case the conjectured bound in [19] is not 
saturated, the perturbative expansion parameters are of the form 

/m \ (4— A) 

S=[j) , so that Pt[6] = (A-4)6, (4.16) 

and (e — 3p)/T 4 ~ 5 2 . This means that the interaction measure (as well as c 2 — 1/3) is not 
a constant at leading order: (3t[{£ — 3p)/T 4 ] ~ 5 2 is of the same order as (e — 3p)/T 4 . This 
in turn implies that the model does not have an effective Chamblin-Reall dual description. 

The above examples suggest that Buchel's bound could be perturbatively saturated when- 
ever (5t[{s — 3p)/T 4 ] is subleading, i.e. that, more generically 

C - ~ 2 (\ ~ c\ U t^lll \ T 9 _ 4 1 ( e _ 3p) , (4.17) 
r] \3 8 J sT [ or J v F> ' y 1 

for some model-dependent dimensionless function f(5, A, T). 12 This proposal, modulo a slight 
modification needed to accommodate the results collected above, is analogous to those pre- 
sented in [44] by means of exact sum rules and a certain assumption for the spectral density 
of the trace of the stress energy tensor. Notice that our formula (4.17) is meant to apply 
to four-dimensional strongly coupled plasmas with massless flavors having a two-derivative 
five-dimensional gravity dual, in a regime where any possible "deconfining temperature" T c 
can be neglected. 13 
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A Details of the calculations 



Tensorial perturbations 

From the action (2.9) it follows that the equation satisfied by the tensorial perturbation 
Zt = 7i. X y is 

^'(^H + J^f)^ ' (A ' 1) 

With the ansatz (3.6) one can check that the only non-zero term in the solution normalized 
to one at the horizon up to first order in \h y d is -2^'° = 1- At second order in \h y d the 
solution is too lengthy to be reported here but straightforward to obtain. 



Scalar perturbations 

The equations for the scalar gauge invariant fluctuations are the relevant ones for the sound 
channel, giving the dispersion relation in (1.2). They are just combinations of the equations 
in [29]: 14 

= 2Hf z OM + A^Hg OM -(l-^^)Hf a OM + 2^H* OM + 
uj tz V u 2 c' x c x y oo 2 c\ u 

+ f^ E + 1 log' £*) H r E ° M + ~H* OM , (A.2) 
V q c l T uj c T J 

n iEOM _ 'P'b ttEOM , LjJ( W c ^X ( \>"b + 0g( c X ~ c xc' x )} ( jtEOM , g C T ttEOm\ 

■ - 9 log' 4 aa c' x (q 2 c T c' T Cx + 2q 2 c 2 T c' x -3uj 2 c x c' x )\ rt uc 2 x rz )' 

In these expressions, represents each of the scalars f,w,$ (the form of their equation is 
the same) and 4> B their background value. Moreover, the notation (j) EOM (and H^ z OM and 
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H aa in [29] corresponds to our H±±. 
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so on) stands for the corresponding equation for the scalar (and the fluctuation H zz and so 
on) in section 3 of [29]. For the coefficient S we have 



_ ^q 3 rj{q 2 -uj 2 ){r 4 -rj) 2 

'ru 2 [q 2 (r A h - 3r 4 ) + 3r 4 cu 2 ] ru 2 \q 2 [r\ - 3r 4 ) + Zr^uj 2 ] 2 ^ ' 



16gr£(2g 2 - 3a; 2 ) 



(A.3) 



We give here the solution to the non-zero fluctuations entering in the sound channel, 
satisfying the normalization at the horizon and Dirichlet conditions at the boundary in the 
case of Z s 



7 o,o _ 
z, s — 



Z 



0,2 



logp 



Z 1 ' 2 = 

v 144^(1 - p 4 ) 



12(1 -p 4 ) ' 

vr 2 (p 4 - 1) + 24(p 4 - 1) log 2 p - 12 logp (4 + (p 4 - 1) log (1 + ip) + 



(p 4 - 1) log + p)(p 2 - 1)]) - 3(p 4 - l)L* 2 (p 4 )" , 



where p = r/r h . We do not report the expressions for all the q 2 coefficients of the solutions 
because of their very lengthy form. 



Vectorial perturbations 

The equation in this channel reads 



IT 



2 C T 



'V 



0. 



For the vectorial fluctuations we can solve order by order with the scaling 15 to 
<l —> ^h y d<\, imposing regularity at the horizon. The result is 



'V 



C 



v 



.■ tPT 



-5+ l-<- 



1-4 (l + e, + e 



+ 0(tt>,q 2 



(A.4) 



(A.5) 



From Dirichlet conditions at the boundary r 
the dispersion relation 

uj = 



00 we can read off the shear viscosity from 



i ^ 2 + 0(g 3 ) . 



(A.6) 



This calculation is summarized in the membrane paradigm formula given in [45] , and it gives 
the well-stated ratio rj/s — l/(47r) with corrections in powers of r^/V* — > 0. 
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Scaling the frequency also with vo — > XhydW gives the same answer, as is the case in [31]. 
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